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F INTRODUCTION 


The field of three-dimensional (3D) underwater acoustics has been 
developed in an attempt to describe the properties and characteristics of 
sound propagation in the heterogeneous ocean environment (Lynch and 
Chiu, 1989). At present, there are three major approaches used in studying 
the 3D underwater sound field. One of these, ray theory, is premised upon a 
high-frequency approximation method. This method gives a geometric- 
optics solution which involves ray-tracing in a spatially varying sound speed 
field. The Hamiltonian Acoustic Ray Tracing Program for the Ocean 
(HARPO) is a versatile numerical code developed by Jones et al. (1986) for the 
computation of 3D rays. The primary deficiency of this theory is that it cannot 
adequately address the behavior of low-frequency sound due to the neglect of 
sound dispersion and diffraction. 

The second approach uses a parabolic approximation to the acoustic wave 
equation, which was introduced by Tappert (1977). A 3D numerical parabolic 
equation (PE) model was developed by Lee et al. (1988) using an implicit finite 
difference scheme. Another 3D PE model was arrived at by Baer (1981) 
utilizing a split-step Fourier algorithm. In a recent study, analytic solutions to 
the 3D parabolic equation were obtained by Seigmann et al., (1990). These 
solutions are valuable for testing the accuracy of 3D numerical models. 

A third approach is the normal mode method. Pierce presented a three- 
dimensional version of this method in 1965. The normal mode method is 
based upon a separable solution to the wave equation. Pierce assumed an 


adiabatic acoustic environment which leads to the neglect of coupling 


presented by Chiu and Ehret (1990). This later 3D normal mode model 
accounts for both horizontal refraction and model coupling. 

Any improvement in acoustic modeling must be quantified. We shall 
endeavor to test the accuracy of the 3-D coupled normal mode model of Chiu 
and Ehret (1990) in regard to horizontal refraction. 

To examine the accuracy of the Chiu-Ehret model, we compare the results 
of their model with those from analytic solutions to the parabolic equation 
arrived at by Seigmann, et al. (1990). The results of the Chiu-Ehret model and 
the analytic solutions of the parabolic equation are compared for two cases of 
horizontal variation in sound speed. As a means of comparing the results, we 
examine the slow variations of the complex pressure envelope function and 
transmission loss. 

Our simulated acoustic field has a range of 100 km with azimuth from 30° 
to 180°. Two variations in the horizontal sound speed field were used. In 
Case I the sound speed varies only with the azimuth angle while Case II varies 
both with the azimuth and radial. 

The normal mode theory used by Chiu and Ehret is described in detail in 
the Appendix. The appendix also includes a description otgiiiegae 
approximation and the analytic solutions developed by Seigmann et al. 


(1990). 


If. ANALYSIS AND METHOD OF MODEL COMPARISON 


The three-dimensional coupled normal mode model was applied to two 
sound speed fields for which exact analytic solutions to the PE approximation 
are available. These analytic PE solutions for vertically invariant sound speed 
were developed by Seigmann et al. (1990). Their development was based on 


expressing the “time-independent” acoustic pressure component as 


fee salen > 
®pr(r,6,z) = U,(7,@)sinl y;z eee ey | . (i -2) 
pe(?, 8,2) = Uplr )sin(y, Jexp 2 = | a | (| exp} 7 


where 


; 
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Co 

ee 
i -|3 AB = 0,1,2,...) 
C, = reference sound speed 


f = acoustic frequency 
H = ocean depth. 


The reduced complex envelope function 
U,,(r,8) = A,(r,@)exp{iO(r, 6)} 


U,, is governed by 
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C(7)= sound speed as a function of position 7. 
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Substitution and separation of real and imaginary parts gives two coupled 


equations governing A, and ©: 
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If one chooses O(7,6) = a(r)@, equation (2.2) may be solved for the amplitude, 


1.€., 





A,(r,0) = F op a | (2.4) 

Choosing a@(r) and the functional dependence F results in a reduced 
envelope function (and ultimately acoustic pressure) while equation (2.3) can 
be solved for the index of refraction n(7) and thus the sound speed field C(7). 
Equation (2.3) is for a wide angle PE approximation. This wide angle form 
will be used throughout. To judge the accuracy of the Chiu-Ehret model in 
calculating horizontal refraction we compare the numerical coupled mode 
results with the analytic PE results for two sound speed fields or two 


functional definitions of a(r). 


A. INDEX OF REFRACTION 


1. Case! 


For Case 1 we choose = Q@(r) = Qk r, or 


where Q is an arbitrary constant. From equation (2.4) the amplitude is 


A,,(1,8) = F(@-a,lnr). 
Let the function F be multiplied by B)ky. The amplitude has the final form 
A,,(7,8) = Bok )(@-aolnr) (2.6) 


where f , is an arbitrary constant with units of length. The square of the index 


of refraction from (2.5) and (2.6) is 
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The resultant sound speed is shown in Figure 1 for a, = 0.005 and 
C, = 1500 m/s. 
2. Case Il 
For Case II we choose 


| 
O = apkor’é (2.8) 


and hence 


A (7,8) = Bok 9(@-Qok gr). Ze) 


The square of the index of refraction is now 
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The resulting sound speed is shown in Figure 2 for @ = 4x10’ and Cy = 1500 


m/s. 


B. ENVELOPE FUNCTION 
We now have two analytic sound speed fields and the associated acoustic 
pressure fields derived from the parabolic equation approximation. 


Considering only j = 1, Seigmann et al., (1990) define pressure as 
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For Case I, using the amplitude and phase expressed in (2.6) and (2.5), 


respectively, we obtain 





Orel, 8.2)=[ 


7 ze 1 D 
xXexp eee | ae (2212) 
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For Case II, using (2.8) and (2.9), we get 





Figure 2. Case II Sound Speed Field for a = 4x10’ and Cy) = 1500 m/s 
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In terms of normal modes, pressure can be expressed as 
C7 
Oyy(r,0,2) = Duplo) nr, 8)ar 7 (2:1, 0) (2.14) 
r 
n 


wilene! 7 eissine n‘” normal mode (eigenmode), k, is the corresponding 
horizontal wavenumber (eigenvalue) and U, is the corresponding slowly 
varying modulation envelope. See the appendix for details. 

Since we are considering depth-invariant sound speed fields, an analytic 


solution is available for the eigenvalues and eigenmodes. The first mode is 


2 
Z,= H sin( 7,2) where Yj = fk? ake a 


Zi 


and thus, keeping only the contribution of the first mode, (2.14) becomes 
Pym (r,@,Z) = (27 U(r, Ae Ho ee sin( 712). (2:15) 


The normal mode model entails a numerical calculation of the slowly 
varying envelope function U,. Therefore, the quantification of model 
accuracy can be achieved by comparing the numerical U, to those derived 
from the analytic PE solutions. The analytic envelope function U,(r,6) for 


Case I is derived by equating (2.12) to (2.15). The resulting expression is 
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With the index of refraction n given by (2.7), the range integral of k, can be 


expressed as 
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N 
1 koro-= | 
U(r,é) = ik 4 (2.17) 
\ ko | 
where 
z 
Z 





: ar 2 i a 4 
g=1-— J eal + Aj99-—43—.,/441+—| — = 400 + 206 + LL 
2\ ko 4\ ko 2\ k ko 


For Case II, the wave number integral, using (2.10), becomes 
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The analytic envelope function U, is then 
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where 


(2.18) 


For Case I we use the analytic envelope function (2.17) to define 


the initial condition at r = ro for the numerical model run. The condition is 


: TT 
1 (komo) 
Ui(10, 8) = 4/7 Boko(9- a In 7 Je ; 


where 
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(2.19) 


For Case II we use the analytic envelope function from (2.18) at r=ry. The 


Case II initial condition is 


: Tv 
H i kor09-2 
Uy (19, 8) = a, Pokol® — aokoro Je ; (2.20) 


where 
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III. MODEL RUNS 


Our computational domain is in a cylindrical coordinate system. The 
horizontal range is from 1 km to 100 km. The azimuth range is 30° to 150°. 
The depth is 4000 m. The domain is the same for both cases. 

We use a frequency of 50 Hz and source depth of 1000 m. The reference 
sound speed is 1500 m/s. We consider mode 1 only. The vertical boundary 
conditions for the problem are pressure release surface and rigid bottom. Two 
range and azimuth dependent sound speed fields as presented in the last 
chapter are used. For Case I sound speed ranges from 1481 m/s to 1496 m/s, 
varying only in azimuth with a gradient of 0.128 to 0.130 m/s/degree (see 
Table 1). For Case II sound speed varies in both range and azimuth (see Table 
1 again). Here sound speed goes from 1437 m/s to 1500 m/s with a radial 
gradient of 0.13 to 0.657 m/s/km and an azimuth gradient of 0.004 to 0.431 
m/s/degree. Both selections of sound speed fields contain realistic gradients 


as observed in the ocean. 


TABLE 1. THE AZIMUTH AND RADIAL SOUND SPEED VARIATION 


Case Sound Speed [nv/s] Azimuth Variation Radial Variation 
aC dC m 
55 F/deg] 5 /km] 


LS Ea 20 50s 
ai 






Range 
r[{km] 










i _[ 1496 _[ 1488] 1481 [0.1303 | oreo fois] 0 | o | 0 

1 [730 [1496 1488_[ 1481 _[ 0.1303 [0.1289 | 0.1275] 0 | 0 | 01 
[100] 1496 [1488 | 1481 [0.1303 | 0.1289 127s 0 | 0 | 0 

1 
[100 | 1487 [1462 | 1437 | 0.4265 | 0.4127 | 0.4001 | 0.1301 | 0.3763 | 0.5992 


WW 


Eigenmodes and eigenvalues are calculated using finite difference 
approximations to (A.8) in the appendix with a 40 m vertical grid spacing and 
a matrix eigenvalue solver. 

To obtain the numerical envelope function a first order differential 
equation ((A.24) in the appendix) is integrated. The radial integration step 


size is 0.5 km. The 61 integration paths were separated by two degrees. 


ilk 


IV; KESUERS 


A. CASEI! 

In comparing the numerical normal mode (NM) and analytic parabolic 
equation (PE) results we focus on the slowly varying modulating pressure 
envelope. The removal of the rapid oscillations, which have wave lengths of 
order 2m/ky, makes the displays of results easier to interpret. 

Figures 3 and 4 show the amplitude and phase, respectively, of the 
envelope function calculated by the normal mode model using the Case I 
sound speed field. Figures 5 and 6 are the amplitude and phase, respectively, 
of the envelope based on the analytic solution to the parabolic equation. 


Figure 7 shows the relative difference, defined as 


Up,( analytic solution) - U,;,, (numerical solution) 








Up, (analytic solution) 


The percent difference is everywhere less than 5%. The difference in NM and 
PE phase is shown in Figure 8. 
Transmission loss at a depth of 1000 m calculated by the NM model, 


Thy wy, is shown in Figure 9. The magnitude of the analytic PE pressure, from 


ee 
[a= fap sin 7)2) (3.1) 


Substituting (2.6) in (3.1), we get 
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Figure 4. Phase of U,,,, for Case I with @ = 0.005 and B = 0.1 
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Figure 6. Phase of U,, for Case I with a@ = 0.005 and B = 0.1 
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Figure 7. Amplitude Difference in Percent for Case I with @ = 0.005 and 
B=0.1 
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Figure 8. Phase Difference in Degree for Case I with @ = 0.005 and B=0.1 
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The analytic PE transmission loss TLp,; is thus 
T Lpg = -20logio|P| 


= 10logi 2" - 201021 0|Boko (8 - ag Inr) - 20log 
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sin——zZ 
Ziel 





Figure 10 displays TLp, at a depth of 1000m. Both transmission loss results 
show that the azimuth variation of sound speed produces lower loss at large 
azimuth angle. 

The difference between the analytic PE transmission loss TLp, and the 
numerical NM transmission loss 1L,),, 1s shown in Figure 11. OGime@eaaine 
modal transmission loss is a function of amplitude alone this difference has 
the same shape as the relative error in the amplitude of the envelope function. 


The difference in transmission loss is everywhere less than 1 dB. 


B. CASE II 

In the second case, we examined a sound speed field that varies in range 
and azimuth, closely representing an eddy or ring structure in the real ocean. 

Figures 12 and 13 show the amplitude and phase, respectively, of the 
envelope function U,,,, from the normal mode model. Note that azimuth 
variation of the amplitude is larger than range variation. Figures 14 and 15 
are the amplitude and phase of the envelope function Up; based on the 
analytic solution to the wide-angle PE. Figure 16 shows the relative 
difference of the amplitude. The percent difference is everywhere less than 
2%. The difference between the NM and PE phases is shown in Figure 17. 


The difference of phase is everywhere less than one degree. 
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Figure 9. TL,,,, at a Depth of 1000 m for Case I with a@ = 0.005 and B= 0.1 





Figure 10. TL), at a Depth of 1000 m for Case I with @ = 0.005 and B= 0.1 
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Figure 11. Difference in TL in dB for Case I with a@ = 0.005 and B= 0.1 
The NM and PE transmission losses, TL,j,, and TLpg, at a depth of 1000 m 
are shown in Figures 18 and 19, respectively. They are nearly the same as can 


be seen in Figure 20 where the error is everywhere less than 0.2 dB. 
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Figure 13. Phase of U,),, for Case II with a= 4x10” and p=0.1 
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Figure 15. Phase of Up, for Case IT with a@ = 4x10” and p= 


a2 
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Figure 17. Phase Difference in Degree for Case II with a= 4x10” and B=0.1 
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Figure 18. TL,,, at a Depth of 1000 m for Case II with @= 4x10” and B=0.1 





Figure 19. TLp, at a depth of 1000 m for Case II with a = 4x10” and p=0a 
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Figure 20. Difference in TL in dB for Case II with a@ = 4x10’ and p=0.1 
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V. DISCUSSION 


In this chapter we address two issues. The first is a quantification of 


horizontal refraction. The second is an assessment of model accuracy. 


A. HORIZONTAL REFRACTION 

Critical to this accuracy test is the use of sound speed fields that cause 
significant horizontal sound refraction. To confirm that the selected sound 
speed fields do cause significant horizontal refraction, we compare the 
corresponding Nx2D NM solutions (the calculation is divided into N vertical 
slices and a 2D NM model is used in each radial direction) with 3D NM 
solutions. 

The envelope functions of the N x 2D solutions have been calculated 
from (A.25) in the appendix. The amplitude of the NM envelope function is 
the same for both the Nx2D and 3D solutions in Case I and Case II. However, 
the phases are different. Figure 21 shows the phase of the Nx2D envelope 
function for Case I. The phase is essentially constant at -45°. In Case II, the 
phase of the envelope function of the Nx2D solution is also constant and its 
value is also -45°. 

Now let us compare the phases calculated by the Nx2D method with the 
phases of the 3D solutions for Case I and Case I, as displayed in Figures 4 and 
13, respectively. Both cases show a 13°-15° phase difference over a range of 
100 km, implying that the azimuthal sound speed variation used is large 


enough to induce sound propagation out of the vertical plane. 
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Figure 21. Phase of U,,,, for Nx2D for Case I with a = 0.005 and B = 0.1 


we ACCURACY 
1. Phase 
In this section we discuss the differences between the analytic PE and 
numerical NM solutions. In Case I, the phase difference between the analytic 
PE and numerical NM solutions is less than 2.8 degree everywhere. The Case 
I] phase differences between the numerical NM and analytic PE solutions, as 
shown in Figure 17, is everywhere less than 1 degree. A quantification of the 


phase differences for the two runs is given in Table 2. 


FABUE 2. PHASE DIEFERENCE 


P| min | imax | min | max | max 
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2. Amplitude 
Results are shown in Table 3. Envelope amplitude of NM and PE 
generally agree well with the largest relative difference near 5%. The largest 
difference occurs in a region that has the largest azimuth gradient in sound 


speed. 


TABLE 3. AMPLITUDE DIFFERENCE 


Difference | Difference [%] 
| min | imax | min | omax | max | max 


3. Transmission Loss (TL) 











The difference in transmission loss between the analytic PE and 
numerical NM solutions is shown in Table 4. The difference for both cases is 


everywhere less than half a decibel. 


TABLE 4. TRANSMISSION LOSS DIFFERENCE 


pT min max min | imax | max 
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VI. CONCLUSIONS 


In this work we have tested a three-dimensional numerical acoustic 
coupled mode model. We have sought to describe its accuracy in regard to 
the modeling of horizontal refraction. In order to test its accuracy, we 
compared the coupled mode model results with two different analytic 
solutions to the parabolic wave equation. 

We have concentrated on the accuracy of the envelope function and 
transmission loss calculation. For the sound speed fields described in Table 
1, we found that the phase error of the slowly varying envelope function is 
lower than 2 degrees and that the envelope amplitude agrees to within 5% or 
better. Our test cases have indicated that the normal mode model agrees 
closely with the analytic solutions. 

We have only used depth-invariant sound speed fields for this test. Asa 
result, the accuracy in modeling mode-mode interactions is not tested here. 
Future tests should be directed at examining the accuracy of calculating 


mode-coupling effects using depth varying sound speed fields. 
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APPENDIX 


A. THREE-DIMENSIONAL COUPLED NORMAL MODE METHODS 
1. Overview 
a. Three-Dimensional (3D) Helmholtz Equation. 
The acoustic pressure P(7, tf) in the ocean is governed by the wave 
equation: 


1 3d7P(7,t) 
C*(F) at? 





V*P(F,t)- =) (A.1) 


where C(r) represents the speed of sound propagation. The cylindrical 
coordinates are 7 = (r, 0,z) where r is range, @ is the azimuthal angle 
measured positive counterclockwise, —z is depth, and V~ is the Laplacian 
operator in cylindrical coordinates, 
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Inclusion of an acoustic source q(7,t), modifies Equation (A.1) to 


Dede 
V2 p(t) =p amt) (7Es) 


For an acoustic point source located at rg with time-harmonic circular 


frequency w, Equation (A.3) becomes 


2 Ae 
Vn(7,t) - = = —476(7 - omen (A.4) 


If we let 
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p(7,t) = &(r)e (A.5) 


where (7) is time-independent, then a Helmholtz equation results: 


V* (7) +k? (7) @(7) = —426(7 -7,) (A.6) 
Ww 
where the acoustic wavenumber is k(7) = FP) For r# r., the Helmholtz 
Equation 1s 
V? (7) +k* (7) B(F) = 0. (A.6a) 


b. Coupled Normal- Mode Solution 

Adiabatic normal mode theory was applied to a range-dependent 
medium by Pierce [1965] and Milder [1969]. The pressure field was expressed 
as a linear combination of the local modes whose coefficients were obtained 
from a system of decoupled ordinary differential equations. The physics of 
coupled normal mode theory includes non-adiabaticity and results in a 
coupled system of differential equations. The theory used in the 
development of the Chiu-Ehret 3D coupled normal mode model follows. 

The first step in the development is the expansion of the time- 


independent pressure component in terms of the local normal modes, i.e., 
J ] 
in Pal 0)Z,,(z:r, 8) (A.7) 
n 


where 


Z, is the n‘” local normal modes at point (7,6). 
ieeis the n'” mode amplitude function at point (1,6). 
c. The Local Normal Modes Z,, 


The local normal modes obey the following equation: 


3] 


2 + (K? ~ i) = (A.8) 


where k, is the horizontal wavenumber associated with the n'” mode at each 
horizontal location. 
The idealized upper and lower boundary conditions are 
Ze Ono): (A.9) 


d 


5, 2nlz = —Hir,8) = 0. (A.10) 


Where H is the ocean depth. The local normal modes can be normalized 


according to the orthonormal condition 
fa 
eae tore (A.11) 
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d. The Mode Amplitude Function P,, 
To obtain the governing equation for P_, substitution of (A.7) in 


(A.6a) is necessary. The use of (A.8) and the farfield approximation 
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following the substitution results in the following equation: 
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The next step is to multiply (A.13) by Z, and then integrate over Z, i.e., 
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A subsequent application of the orthonormality condition (A.11) on (A.14) 
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Defining coupling coefficients as 
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the coupled system (A.15) can be recast as 


0 1 3 d Igoe 
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e. The Envelope Functions U,, 
Following the work of Chiu and Ehret (1990), the mode 
amplitude function P, can be separated into a slowly varying envelope 


function U, and a rapidly varying component ae 1:€., 
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P, =U, (r, 0)etonl"9) (A.17) 


On [ kn(r, Oar. (A.18) 


A substitution of (A.17) into (A.16) yields 
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which can be expanded as 
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Dividing (A.20) by i2k e °n gives the governing equation for the envelope: 
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and ? and @ are the unit directional vectors in r and @, respectively. 
f. Transmission Loss TL 
The mode amplitude function P, is related to the envelope 
function U, and the phase @, by (A.17). The local normal modes Z, are 
solutions to the eigenvalue-eigenfunction problem (A.8)—-(A.11). Once P,, and 
Z,, are computed, pressure can be calculated from the product of P, and Z,. 


At each point (r, @, 2), p is 
©? (r,6,z) =[Re(®)]° +[Im(@)] (A.22) 


where 


Re| ®(r, 6,z)| = LIRe(Frn a (z) 


and 


m(®) = ¥'[Im(Py )]- Zm(2) 


In (A.22), Re and Im are used to denote the real and imaginary 


parts, respectively. Transmission loss can be computed as 


TL = -10 log,)®(r,6,2). (A.23) 
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2. Numerical Solution 
a. Numerical Model 
To obtain a normal mode solution we must solve the equations 
for the eigenvalues and eigenmodes (A.8)—(A.11) and the equation for the 
modulation envelope (A.21). These equations can be solved numerically. 
b. The Eigenvalues k, and the Eigenmodes Z 
Equation (A.8) is approximated using central finite differencing. 
This approximation casts (A.8)-(A.10) into a matrix algebra problem. The 
eigenvalues and the eigenmodes of which can be determined using the 
iterative OR method (Acton, 1970). 
c. The Iterative Method 
To obtain the envelope function U,, we need to solve equation 
(A.21) which can be rearranged to form a set of first order partial differential 
equations (PDE) with smaller terms put on the right-hand side. We can solve 
for U, in an iterative fashion. The iterative equation is 
2 2 


NT 


Oy ere ele alee 


Hi 


where ie is the solution at the i” iteration for mode n. This set of first order 
differential equations is integrated using a Runge-Kutta method of order 5 
and 6 (Acton, 1970). 
3. The N x2D Method 
In the N x 2D method, the sound channel is divided into N vertical 


slices and the two-dimensional normal mode solution is solved for each slice. 
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It combines the results in each vertical plane to construct a 3D field. This 

method provides an approximation to a 3D solution for a 3D sound field. 
Diiicmim eden slice the envelope function U, 1s independent of 

azimuthal variation effects. With the azimuthal terms ignored, the governing 


equation becomes 


ul! -] 
0 ; ; ON 2 i] 
2 : EY fut ; y AU, 7 Dy wa ss ys Cele (A.25) 
: M1 Or 1 
where By — 
i 
1 ok 
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a 1 x = 
Can ™ 5 Yn ne One 
ae 
or 
1 7 (dm - On) 
: ap anak si Bee (m # 71) 
Ayn = j : 
Smet (m= 1) 
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B. THREE-DIMENSIONAL PARABOLIC WAVE EQUATION 
1. Parabolic Wave Equation (PE) 
The following work was developed by Tappert (1977) and Lee (1988). 


In cylindrical coordinates, the Helmholtz Equation (A.6a) becomes 


eb 100 100 Oo 


me i 








where 


ky = W/Co 


o7 


C, is a reference sound speed 


w is the harmonic source frequency 


n is the index Of retraction 


n=n(r,0,z)=Cy /C(r,@, 2). 
The PE approximation assumes 


@D(r,6,z) = u(r,8,z) Ur) (A.26b) 


Where v(r) is rapidly varying and u(r, 8, z) is a modulation. 


By substituting (A.26b) into (A.26a), one obtains 


671 }3s 200) 6u 1 Ofu Ot Or 0 Von 
= +t et kon’ ujo+| —> +—-—|u =0. (A.26c 
a. r vdrjer r*@0* Gz? : | Ss r a | | 


If v is determined by 


Z 
—+-—+t+kjv=0 A-2Z7 
a ( ) 


then u must satisfy 


en (1 26v)eu 102" O2u 2, 2 
—— +| —+ —— |— + ~—~ + —— thin” - Iu = 0. A.28 
ar? 2 v =o re ao? Oz a | ( ) 
The solution to (A.27) is 
ple 
u(r) = H (kor) ca 2) expli kyr “|| (A.29) 
tkor 4 


Ee : 
where Hog is a Hankel function of zeroth order. 


By substituting (A.29) into (A.28), one gets 
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Q2 ae | ele 
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Roe Oa Seidl = 0 (A.30) 
y 


which can be factored into 
QO. . Ca . 
— t+iky —ikgQ | — + tk + tkypQ ju =0 rol) 
Or Or 


where 


1 

Z 2. 

9 lec 1 da {2 

See eal ee res el 
. | | | ké dz? kér? = 


Q is called the square root operator. If one considers only outgoing waves, the 


equation to be solved is 


(= + ikg - ky — (6) (A.32) 
Or 


Let 
| 
oO (1 +X + Y|2 (Ac33) 


SO tnat 

aloe 

ké az? 

melon 
ker? 067 


X =(n?-1}+ 


and then expand Q in a Taylor series: 


1 
Q=[1+X+Y]2 = 14+5X-2 x7 4 Y4... (A.34) 
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Inclusion of the first four terms in the expansion leads to a “wide angle” 


version of the PE approximation. Substituting (A.34) into (A.32) gives 


2 
Ou iky |, 2 107 al 1 0 less 
— = — [n° - 1) + —— - — h(n -1)+—S— st a eee. (A. 
ore eed kg Oz" | | a k6r? 607 = 


2. The Envelope Function U,, of the Parabolic Equation 


Here we present the analytical results of Seigmann et al., (1990). Their 
development uses the pressure release surface and rigid bottom boundary 


conditions, 1.e., 


Se ied ee (A.36) 


A solution to equation (A.35) for a depth-independent sound speed field and 


subject to the boundary conditions stated above can take the following form: 


z Z 
mie ‘ 2 ik gr Yi 1 ial 
Wi to 2). = U,(r,8)sin(y ;z)exp Aa a) c ar (A.37) 


where 


-[i3]4 


H is the ocean depth, and U,, is the envelope function of the parabolic 
solution. 


By substituting (A.37) into (A.35), one obtains the equation governing 
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(A.38) 
Letting 
Uy (7,0) = Ay(r,aye 7) (A.39) 
where both A (7,8) and O(r,@) are real quantities, (A.38) becomes 
2 : 
0A IkjA ly 1 Dailies 
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A separation of the real and imaginary parts of (A.40) yields 
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Given A, and ©, a pressure field is determined, and the index of refraction can 


be calculated from (A.42). The resultant sound speed field can then be used 


in the normal mode model. 
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